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ABSTRACT 


Selective problem areas relating to computer based ESM 
system software development are discussed. Digital simula- 
tion is presented aS a practical approach to solving software 
development problems. An algorithm to simulate a hard-wired 
EW signal preprocessor is developed. An application oriented 
high level language version of the algorithm is presented. 

A FORTRAN program which generates realistic data for 
Signal analysis is described. Using the simulated data, an 
instruction level simulation program, and @ language compiler, 


the algorithm is tested using an IBM 360 computer. 
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I, INTRODUCTION 


A. STATEMENT OF THE PROBLEM 

Electronic Warfare Support: Measures (ESM) is that 
Givision of Electronic Warfare (EW) concerned with exploita- 
tion of enemy use of the electromagnetic spectrum. ESM 
includes actions taken to search for, intercept, locate, and 
analyze electromagnetic radiations to support tactical 
military operations. 

The continuously increasing density and complexity of 
today's electromagnetic environment, created by friendly as 
well as hostile forces, has fostered the requirement for more 
effective performance and greater sophistication in operational 
ESM systems. Current receiver technology has kept pace with 
this escalation, and signal intercept is virtually assured. 
However, the resulting increaseca data rates, existing and 
forecast manpower limitations, and limited reaction time 
available to the operational commander have imposed a require- 
ment for faster and more accurate Signal analysis, classifica- 
Pilon, dna aitexrpretation. Currently, major efforts are being 
put forth to achieve an effective ESM system design incoxrpo- 
rating high speed digital computers for the purposes of 
automatic signal analysis and system control. 

To perform effectively in present and future environments, 
an automated ESM system should exhibit the following 


GijvabacteLiStLes: 





Me selec mrrcauency range . 


Simultaneous coverage of entire frequency range. 


Dame rOoapLLity Of intercept . 


Low false alarm rate, 


Capability to recognize pulse or radar-~like signals. 


NAN MO S&S W NH 


Capability to determine parameters associated with 
radar pulse trains. 


7. Capability to classify signals according to type and 
source. 


Mo ween tom rinaingd Capability. 
9. Capability to handle data rates. 


10. Capability to detect signals at ranges to radar horizon. 
The successful design and development of a system exhibiting 
these characteristics rest upon three equally important 
elements: system hardware, system software, and the man- 
machine interface. 


l. ESM HaxrdwareysSystem Description 


8 eS ee a Re RS RE 








A @iscussion of ESM hardware as a specific problem 
area 1S beyond the scope of the thesis. A system concept 
compatible with the state of the art approach 1s presented 
here aS a basic foundation for this study and future efforts 
in the ESM arena. 

Conceptually, an airborne ESM system with the charac- 
teristics previously described could be implemented utilizing 
a combination cf Instantaneous Frequency Measurement (IFM) 
and Superhet receivers. The IFM’receivers, operating in 
conjunction with a high-gain pencil beam antenna, would cover 
the frequency bands of interest and function as the primary 
fe@Ulare@lenmmiaeced mrs On the System. In addition to contrib- 


uting to increased system sensitivity, the high-gain pencil 





beam antenna functions to reduce the inStantaneous data rate 
aeeyeemecOngunction With azimuth centroiding, provides an 
accurate means of determining Signal direction of arrival. 

Additional data rate reduction is accomplished by 
routing each IFM output through a Pulse Train Separator (PTS), 
which functions as a signal preprocessor to separate stable 
pulse trains from the received data. Parameter measurements 
are performed on those stable trains separated and a computer 
word is generated to tag the associated emitter. The required 
capacity of the PTS in number cf pulse trains separated is 
considered paramount to system operational effectiveness. 

A computer word is generated for each pulse associated 
with non-stable pulse trains and those stable puise trains 
exceeding the capacity of the PTS. These data are routed 
through a residual channel of the PTS to a high-speed digital 
computer for deinterleaving and parameter measurement. The 
parameters associated with each signal are then passed to a4 
Main frame computer for classification, threat analysis, and 
tactical environment update. 

The IFM system employs an omni antenna for side lobe 
Suppression to prevent Signal poke~-through in the receiving 
antenna minor lobes. The IFM sensitivity and antenna gain 
provide a capability for radar detection cut to the radar 
horizon for aircraft altitudes in the neighborhood of 
30,000 feet. 

Priority Signals, unclassifiable signals, and signals 


too complex for automatic enalysis are automaticaily assigned 
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to digitally controlled superhet receivers for refined 
analysis. The superhet receivers cover the frequency bands 
of interest and operate in conjunction with an omni or low- 
gain broad beam antenna. 

2. Man-Machine Interface 

The human operator element of an automated ESM system, 
though normally an area addressed by human factors engineers, 
impacts heavily on the software requirements and will be 
discussed here in that regard. 

Automation does not eliminate the requirement for 
human operators, but it does significantly affect the role 
the cperator piays in the ESM system and it may seriously 
retard the operator's proficiency and ability to function 
efficiently as a back up in the event of a system failure. 
Thus, in @deciding which system functions should be automated 
and which should be assigned to the human operator, it is 
imperative that the operatcr be assigned those functions that 
maximize the use of his capability to exercise judgment, 
improvise, reason inductively, and profit from experience. 

The human operator must be able to control and direct 
the operation of the system. To do so, he must be provided 
relevant information by the system at the proper time, in the 
proper format, and on the proper displays. To obtain efficient 
and reliable software, these requirements must be clearly and 
accurately defined in the planning stage of the ESM system 


development. 


te D 





SB. “Hom S@eftware 

It is felt that software represents the most trouble- 
some problem in the deSign and development of an automatic 
ESM system. AsSSuming cost and risk vary Girectly with 
complexity, ESM software will account for a substantial 
percentage of total system costs and be the riskiest part of 
the system development. Since system performance hinges on 
software performance, implementation of an effective software 
test plan would serve a threefold purpose of cost reduction, 
risk reduction, and increased system performance. 

The ability to successfully implement such a4 test 
plan requires a Cdetailed description of the systems functions 
and a software modular building process. Each software 
module or set of modules is then integrated into the system 
and tested against sets or subsets of system functions 


meeceoraing £O prescribed test plans. 


EpeeeocOrs OF Til THESIS 

Simulation, via computer software, 1s one technique that 
is useful in software development and verification cr testing. 
Realistic data, created by successful simulation of the ESM 
system hardware operating in a dynamic environment, utilized 
as inputs to a high-speed digital computer would provide the 
framework for implementation, testing, and optimization of 
ESM signal processing alcorithms. The inherent ability to 
GOntvolvtie density of the electromagnetic environment through 


Simulation would provide the capability to analyze worst case 
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situations, thereby enabling system designers to estimate 
computer capacity and reserve capacity requirements as well 
as data rate reduction requirements. 

Given these simulation capabilities, an instruction-level 
Simulation of the target computer provides a powerful means 
of dehugging software and aigorithm testing for accuracy and 
speed of overation. Applied Technology Division of Itek 
eempOraiton, sunnyvale, California, has provided the Naval 
Postgraduate School with the programming support system for 
Eneie ATAC computer, Instruction level simulation is only 
a part of this software package which includes the compiler 
Preogtams £om Programming Language for ATAC (PL/ATAC), a high 
level languaye developed for ATAC applications programming. 
Acquisition of the ATAC Programming Support System (APSS) 
gives impetus to the utilization of simulation as a technique 
for ESM software development and testing at the Naval Post- 
graduate School. 

Total hardware simulation of the ESM system previously 
described is beyond the scope of the thesis. Total system 
Simulation does, however, lend itself to a software "build" 
process and as a first effort of the process, three specific 
objectives of this work are as follows: 

1. Provide a digital simulation of an IFM receiver operating 


in conjunction with a high-gain, pencil beam DF antenna 
in a dynamic radar environment. 


2. Provide a digital simulation of a hard-wiredpulse train 
deinterleaver fed by an IFM receiver. 


3. Describe end demonstrate the capabilities of APSS, 
Paiariecwlarly, the LL/ATAC programming language. 





The overall objective of this work is to provide a building 
block as a basis for follow on study of ESM system scftware 


at the Naval Postgraduate School. 
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Ti. COMPUTER SYSTEM FOR ESM SOFTWARE DEVELOPMENT 


eee tee eee ee Ee ee eee eee ee ce mee 


A. GENERAL DESCRIPTION 

APSS was written for Applied Technology's Airborne 
Computer (ATAC). ATAC, specifically designed for Electronic 
Warfare applications, is well suited for ESM software design, 
development and testing. A number cf features of this 16-bit, 
fixed point, two's complement computer permit the establishment 
of a meaningful baseline for determining the relative merit cof 
algorithms applicable to computer based ESM systems. Special 
features of interest are: 


1. 437 nanosecond major cycle time. 


2. Sixteen general registers which may be used as 
accumulators or as index registers. 


3. Microprogremnable processor, permitting optimization of 
the instruction set for specific applications. 


4, Memory types can be mixed in any combination of low power 
Read Only Memory (ROM), fast bipolar Random Access 
Memory (RAM), and slower core memory devices. 


5. Memory expandable to 65K words. 


6. Four basic instruction addxessing modes: 
register, immediate, direct, and direct indexed. 


7. Three supplemental addressing modes: 
register indexed, immediate short, and direct register. 


8. Provisions for up to 1é Direct Memory Access (DMA) 
channels. 


S. Provisions for up to four simultaneously active DMA 
channels. 


10. Twelve hardware priority interrupts. 


it eeiecguliar, compact; construction. 
Additional detailed information, including the instruction 
set, may be ohtained from ATAC, Applied Technology Airborne 


Cormeucen, vol. L, Principles of Opveration. 
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B. ATAC PROGRAMMING SUPPORT SYSTEM (APSS) 

This software support system allows for the preparation 
and debugging of programs via an IBM 360/370 computer system. 
Peeo, through simulation of the ATAC computer and associated 
peripheral devices, permits the applications programmer to 
optimize programs for executicn speed and memory utilization. 
The APSS system is written in a combination of PL/1, FORTRAN 
IV, and IBM assembler language and includes the programs 
listed below. 


APSS Monitor 

ATAC Instruction Assembler 

ATAC Instruction-Level Simulator 
Relocating and Linking Loader 
Microcode Assembler 
Microcode-Level Simulator 

EE ATAC Compiler 

Library Control System 


Subroutine Library 
Arrangements were made with Applied Technology to use portions 
of APSS and consequently, the ATAC Instruction-Level Simulation 
Program and the PL/ATAC Compiler are currently operational on 
the IBM 360/67 at the Naval Postgraduate School. Information 
concerning the implementation of these programs and instruc- 
tions pertaining to the use of the programs are contained 
in Appendix A. Additional APSS information is available in 
ATAC, Applied Technology Airborne Computer, Volume 2, ATAC 
Programming Support System (APSS). The documents in this 
volume provide a detailed description of how to set up and use 


the APSS syStem. Those documents applicable to this work are: 
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Pom aetoeneOL instruction execution. 


2. Tracing option of the 16 general registers. The contents 
Of each register can be printed after the execution of 
cinmnis EeCt 1 Ol. 


3. Accumulated program execution time printed with each 
line of trace. 


4, Instruction execution histograms by frequency of 
execution or time of execution can be created and printed. 


5. Partiai and total simulated memory dump options. 

6. Initialization of program data arrays allowing data to 
be read into virtual program memory prior to program 
Simulation. 

The APSS Instruction-Level simulator program includes a large 
repertoire of error messaces to assist the applications 


programmer in debugging job management and ATAC instruction 


usage errors. 


C. PL/ATAC 

PL/ATAC, Programming Language for ATAC, is a high level, 
block structured language designed specifically for the ATAC 
computer. The programming structure and numerous capabilities 
cf the lanaquage are demonstrated in Appendix C via the Pulse 
Train Separator simulation program written in PL/ATAC. in 
addition to hardware simulation, a specific objective in 
writing this program was to provide programming examples of 
the PL/ATAC languace to assist the interested programmer in 
learning PL/ATAC. The PL/ATAC manual (Volume 2, Section II, 
ATAC) is an excellent reference but may present some difficulty 
if used as a learning text. 

PL/ATAC is a powerful language providing the programmer 


with programming control appxoaching that of assembly language 
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without compromising the high level of the compiler language. 
Additionally, ATAC assembly lancuayge may be intermixed with 
PL/ATAC increaSing programming flexibility and control with 
no serious degradation of program readability. 

An analysis to determine the effectiveness of applications 
programming in the PL/ATAC language versus ATAC assembly 
language was conducted by the author utilizing an assembler 
optimized ATAC benchmark program. The results of this analysis 
indicated a memory utilization increase by a factor of 
approximately 2.3 times the ATAC assembly language requirement 
and an execution time increase by approximately 2.4 times the 
ATAC assembiy language execution time. 

The PL/ATAC programmer loses some fiexibility, as compared 
0 the assembly language programmer, due to the limitation in 
the number of registers available (a maximum of five) to him. 
Also, the PL/ATAC compiler currently installed at the Naval 
Postgraduate School does not have a double precision capability. 
Double precision calculations, if required, must be accomplished 
in assembly language resulting in increased core space and 


execution time required. 
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TIZI. ESM ENVIRONMENT/IFM Sie On 


Investigation revealed the existence of a computer program 
(Ref. 1] designed to Simulate an IFM receiver operating in 
conjunction with an omni antenna in a dynamic redar environ- 
ment. This program was modified to incorporate a rotating 
DF antenna system. Routines were added to generate unstable 
pulse trains in the form of pulse 'jitter' and/or 'sStagger'. 
Documentation was added to facilitate program readability. 

The DATA Simulation Program generates Simulated puised 
radar Signal level data arriving at an IFM receiver via a 
rotating antenna system. Current capacity allows for simula 
mom Of Wo £O ten emitter platforms with up to five ecinitters 
per platform. This capacity may be increased by increasing 
the array dimensions of the parameters associatec with each 
ematter—-platform combination. Emitter platform and IFM 
PreefOorm initial position information in X-Y grid coordinates, 
course in degrees clockwise from North, and speed in knots 
are set by the program input data. 

“jekad (uber Sheaterersilseres, set by the program input date, 
consist of transmitter frequency in MEz, pulse repetition 
materval (PRI) in mitliseconds, pulse width (PW) in micro- 
seconds, transmitter power in kilowatts, antenna gain in db, 
antenna keamwidth in degrees, antenna Scan rate in seconds 
per revolution, and pulse jitter and/or stagger characteristics. 


A sin(x)/x antenna radiation pattern is assumed for ail emitters. 
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IFM receiver characteristics, set by the program input 
data, consist of receiver threshold sensitivity in dbm, upper 
and lower frequency limits of the passhand in MHz, minimum 
time measurement capability in microseconds, antenna gain in 
db, antenna scan rate in seconds per revolution, and antenna 
beamwidth in degrees. 

Additional program inputs include simulation start and 
stop times in seconds, frequency bin size of the IFM digital 
processor, and random number generator calling arguments. 
Complete input data description is depicted in the program 
listing following the description of the program variables. 

A complete program listing is cecntalined in Appendix E. 

Utilizing the input data described, the Data Simulation 
Program computes the posSiticn, course, speed, and range of 
each emitter platform relative to the IFM platform. The 
closest point of approach (CPA) and distance to the CPA is 
computed for each emitter platform to determine if and when 
a platform will be within the range of detection of the IFM. 
Those emitters out of the detection range during the simulation 
are listed as such. The IFM antenna's and all emitter 
antennas' intial pointing directions are set by a random 
process utilizing a uniform random number generator. The time 
at which each emitter platform enters the capture area of the 
IFM antenna's main beam is calculated and used as a basis for 
computing a random transmit time for the emitters associated 


ween that platform. 
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Emitter maximum effective radiated power (P.) is then 


computed utilizing the mathematical model given by 


a 2 2 
on or 


where 
Py = Transmitted Power, 
G, = Emitter Antenna Gain, 


and A = Wavelength of Signals. 


Emitter maximum range of detection (Ro) 1s computed by: 


L. 
a 2 
Rin (PE G//T})) 


where 
e = IFM Antenna Gain, 
and ee LE Receiver Threshold Sensitivity. 


For each emitter within detection range, the minimum antenna 
gaan (S ) required to break the IFM- threshold sensitivity 
is then computed hy 

2 


5 


> min Ty /Gy Pr - 


where 


mo Proeance tosethesemitter, 


eu 





Using the value of emitter transmit time previously calculated, 
the emitter antenna pointing direction is calculated. Emitter 
antenna gain as a function of the relative angle between the 
emitter antenna and the IFM antenna is then calculated and 
compared with the minimum antenna gain required to break the 
IFM threshold sensitivity. If the Signal level is sufficient 
to break the IFM threshold, the time of arrival (TOA) of the 
pulse at the IFM receiver is calculated. 

If the pulse does not break the IFM threshold, the emitter 
transmit time iS incremented by the emitter PRI, the antenna 
position is updated, and the antenna gain calculation process 
is repéated until a Signal level of sufficient strength 1s 
transmitted. Should no pulse be received in the scan time 
of one beamwidth of the fFM antenna, the subsequent scan of 
the receiver antenna is checked and the process continues 
until a pulse is received or until simulation stop time. 

The time of arrival of the first pulse from each emitter 
1S stored in an array and, upon completion of this process 
for all emitters, the array is searched for the lowest time 
of arrival. Tne receiver antenna bearing at this TOA is 
Salctusatead ana provided as output data along with the parameters 
Of the associated emitter. The emitter platform position is 
updated and the time of arrival of the next pulse from this 
emitter is calculated and stored in the array. The array 1s 
tnen searched for the next lowest time of arrival and the 


process repeated until simulation stop time. 
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Miso evewOL et neshata Simulation Program consists of 


parameters associated with each pulse received by the IFM 


receiver in the proper time sequence. These parameters 


include the frequency, pulse width, bearing and time of 


Srrival. 


This data, as prepared for the ATAC simulation programs, 


consists of five integers for each pulse. These integers 


fmemreesent the following: 


ik 


Emitter bearing. An integer in the range of 0-359. 
This word represents the IFM antenna bearing at the 
time of arrival of the pulse. 

Frequency bin. The bins are numbered in ascending 
order of frequency from 1-128. This range may be 
altered as desired by the programmer. 

Pulse width. The pulse width word is an integer 

in the range 0-65.535. The time scaling factor is 
arbitrary. Current convention is that the least 
significant bit of the pulse width word corresponds 
to .l1 microseconds. This allows for representation 
of pulse widths from 0-6.5 milliseconds to within 
Wom mLecOseconds:, 

Time of arrival. Time of arrival is expressed as 

a double precision (32 bit) integer. Word 4 is 

the least significant half of this integer and holds 
time values up to 0.84 seconds with the least signi- 


Pea Vem in ecOrsesPongding to 25.6 microseconds. 





eee wiceortea-rivyain Fhe Woper half of the time word 
Motmeseteime values wp to 15.3 hours. 
The Data Simulation Program produces these outputs in 
printed format and in a punched card format suitable for 
acceptance by the ATAC simulation program. A sample of 
the program output iS contained in Appendix B following the 
meogcam listing. The program flow chart is illustrated in 


Mmroures 1 and la. 
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| INITIALIZE CONSTANTS AND SCALE FACTORS 


— 


READ CALLING ARGUMENTS FOR RANDOM NUMBER GENERATORS, 
SIMULATION START AND STOP TIME, IFM CHARACTERISTICS, 
AND IFM PLATFORM COURSE, SPEED, AND POSITION, 


| INITIALIZE IFM ANTENNA POINTING DIRECTION | 
“OUTPUT IFM CHARACTERISTICS | 








eee READ EMITTER PLATFORM INFORMATION| INFORMATION Bye 


a pe 
< PLATFORM ALT (TUDE 2) Ze En0 > genes 
u : oe 


rh 


OUTPUT EMITTER PLATFORM INFC INFORMATION 


VC gueree se ee Te OE 1 ae ke 7 Se eS re Ue EROCELEA ST 8} sea Ts Ser bee Ce | 


COMPUTE EMITTER PLATFORM COURSE , SPEED, | 
AND RANGE RELATIVE TO IFM PLATFORM 











Mg Vn oo 
ee ee 
COMPUTE BEARING, CPA, AND DISTANCE TO 
CFA OF EMITTER PLATFORM 
ee 
COMPUTE TIME WINDOW IN WHICH THIS PLATFORM 
WILL BE IN RECEIVER ANTENNA BEAMHIDTH 


ee aS foe aan puree eae aT meets epieareee ean semmpeatemeeee ALS Oe a a oe ee ee ae me ee Ss 
READ EMITTER CHARACTERISTICS, CONVERT TO INTERNAL 
FORM, AND COMPUTE MAXIMUM EFFECTIVE RADIATED POWER 


cele gui ya DN, ME ae ee 
| OUTPUT EMITTER CHARACTERISTICS | 


as Siem 


a 
COMPUTE MAXIMUM RANGE AT WHICH IFM CAN DETECT EMITTER | 


—— 

















ee 8 Beas 
| COMPUTE TIME OF ARRIVAL OF FIRST PULSE FROM EMITTER | 
Yes 7. aay 











a = 
IS THIS THE LAST EMITTER ON THIS PLATFORM 


FIGURE 1. Data Simulation Program (Simplified Flowchart) 


NO) 
cr 





gel BEE Uo ls ie 


FIND LOWEST TOA AND OUTPUT PULSE DATA FOR THIS E 





COMPUTE TOA OF NEXT PULSE FROM THIS EMITTER 
AND UPDATE EMITTER PLATFORM POSITION 






SIMULATION STOP TIME 


Yes 


STOP 
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HITTER | 


Data Simulation Program (Simplified Flowchart) 





IV. THE PULSE TRAIN SEPARATOR 





Pa CONCEPT OF OPERATION 

Data received by the IFM are routed to ieee Ewalecie Me Slike: 
Separator (PTS). If a pulse train has a stable pulse repe- 
Beetom trequency {(PRF) the PTS locks on the train and tracks 
the train for a period of time determined by the receiver 
antenna scan rate and beamwidth. The signal's pulse width 
(PS) and pulSe repetition interval (PRI) are meaSured and 
digitized along with azimuth, time first seen and frequency. 
These five parameters are used to tag the associated emitter 
and are routed to the main frame computer for emitter classi- 
fication. Direction of arrival is computed by azimuth cen- 
troiding based on the receiver antenna bearing at the time 
Of arrival of the first and last pulse of the signal. The 
Pe-sapalit? Exists that, due to sin(x)/x characteristics 
of emitter antennas, centroiding will be performed on ragged 
Slonal inputs. This is not considered degrading to system 
performance as meaningful direction of arrival information 
meee Stalie result. 

An inherent characteristic cf the Pulse Train Separator 
is that those steble trains with the lower values of PRI will 
be separated first. This feature of the PTS is also desirable 
as the data rate reduction role of the device is enhanced. 
Paewsegurea capacity Cr the device in this role of data rate 


TeauctlLOoneias not been eStablished and 1S conSidered a most 
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tiipgebeaitenifaageor In the overall ESM system design. Failure 

to provide the proper level of preprocessing could result 

in system software and hardware overloading and, consequently, 
degraded system performance. 

Those stable signals exceeding the capacity of the Pulse 
Train Separator and non=-stable signals such as frequency 
agile signals, jittered pulse trains, staggered pulse trains, 
and other exotic signals are routed through a residual signal 
channel of the pulse train separator to the high speed mini 
computer for parameter measurement. For each of these pulses 
received, a computer word will be generated for frequency, 
meme Or arrival, pulse width, and direction of arrival. Pulse 


einterleaving and such parameter meaSurements as pulSe 


Ox 


repetition interval and pulse jitter and stagger character- 
istics will be performed by the mini computer for emitter 


moses) fication. 


Eee fo OLMULATION PROGRAM 

The Pulse Train Separator Simulation Program is designed 
to Simulate a hardware preprocesSing system with a capability 
of separating we to a maximum of ten stable pulse trains from 
the output of an IFM receiver. 

Input data to the program consist of a series of five 
integer words for each pulse generated by the Data Simulation 
Program; direction of arrival in degrees clockwise from North, 
frequency by bin number, double precision time of arrival 


(two words) in seconds withthe least significant bit 


ZG 
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corresponding to 25.6 microseconds and pulse width in seconds 
Peenethe least Significant bit corresponding to .1 microseconds. 
These data are loaded into a large input huffer prior to the 
eearceot the program Simulation. 

Each pulse within the buffer is compared in frequency 
with previously received pulses within the same receiver 
antenna beamwidth. If no match in frequency occurs the pulse 
is considered a new Signal and a pointer address is stored 
in a New Signal Pointer File. If a frequency match occurs, 
frespulse width of the signals are compared. Upon obtaining 
a match in pulse width the matched pulse's pointer address 
is stored in a five element array within the New Signal 
Pointer File. Upon receipt of five pulses of a pulse train, 
the program calls a local re-entrant procedure (PRICOMP) 
wherein the stability of the pulse train is determined. For 
purposes of this simulation a pulse train 1s considered 
unstable if the variation in time interval between pulses 
is greater than 51.2 microseconds and/or three or more 
missing pulses are encountered in ragged Signals. Upon 
Pec eeR ning a pulse train to be Stable, the pulse repetition 
interval of the signal is calculated and the signal is 
flagged as a stable train. The number of available channels 
in the Pulse Train Separator is then reduced by one and the 
Signal is zeroed in the input buffer as each additional pulse 
is received. Parameters associated with this pulse train 
are stored and the signal bearing is updated as each pulse 


is received through incrementation of the input buffer. 
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Upon receipt of the last pulse expected from the pulse 
imei, the centroid direction Of arrival is calculated and 
the signal information is zeroed in the New Signal Pointer 
File as well as in the first five input buffer lcoations. 
The time of arrival of the last pulse expected is computed 
utilizing the scan time for one receiver antenna bandwidth. 

Non-stable Signals are flagged as such and remain in the 
input buffer. Upon completion of incrementation through the 
input buffer, all stable pulses, except those from trains 
that remain due to PTS capacity having heen exceeded, have 
been zeroed from the input buffer. Parameters associated 
with these trains are stored in sequence according tc pulse 
repetition interval. 

The pulse information remaining in the input buffer, 
associated with the unstable pulse trains and with pulse 
trains exceeding the PTS capacity, is then compressed within 
Mieminiput butter, thus allowing for simulation of the mini 
computer Signal processing algorithm. In the initial design 
concert of this Simulation program, these remaining data, 
representing the Pulse Train Separator reSidual channel, 
would be processed by an ATAC assembly language program 
incorporated as a local procedure within the PTS program 
uSing the code selection head as described in the PL/ATAC 
manual. 

A complete PL/ATAC program listing of the Pulse Train 


separator 1S contained in Appendix C. An assembly language 
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program generated by the PL/ATAC compiler is contained in 
Appendix D. A program flowchart is illustrated in Figures 


Zand 2a. 


Ob 
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FIGURE 2a. The Pulse Train Separator Program Flowchart 
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V. DISCUSSION OF SIMULATION RESULTS 


Pee SUMMARY OF RESULTS 

The Pulse Train Separator as described in Chapter IV 
was Simulated utilizing the ATAC Programming Support System 
and the PL/ATAC Programming Language. The particular hardware 
System Simulated was arbitrarily given the capability of 
deinterleaving and preprocessing up to nine radar pulse trains. 

inoue adata, provytcea by the Data Simulation Program, 
consisted of a mixture of stable and non-stable pulse trains 
from eighteen emitters located on six platforms geographically 
positioned so as to be intercepted in one receiver antenna 
beamwidth. The receiver antenna beamwidth was arbitrarily 
set at five degrees, 

Beeiemcescece! Ver antenna scan rate Of two seconds per 
revolution, the Data Simulation Program generated a total of 
365 pulses as IFM output during one antenna main beam sweep 
of the eighteen emitters. These pulses represent approxi-~ 
mately 28 milliseconds of real time data. The PulSe Train 
Separator Simulation Program, in processing nine of the pulse 
trains, removed 277 pulses from the input data and routed 
88 pulses through the residual channel. These 88 pulses 
represent the unstable pulse trains and the lower pulse 
repetition frequency (PRF) signals. Accurate pulse repetition 
interval calculations as well as accurate azimuth centroiding 


was performed on all preprocessed pulse trains. 
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B. CONCLUSIONS 

The Data Simulation Program provides realiStic data to 
the APSS programs. The program produces reascnably complex 
data in the form of jittered pulSe trains and staggered 
pulse trains. Pulse trains transmitted at sufficiently low 
power levels are produced with missing pulses due to sin(x) /x 
effects. 

The Pulse Train Separator Simulation Program demonstrated 
the use of the PL/ATAC language and functions as designed to 
preprocess ESM signals. This program is somewhat inefficient 
in its present form for two reasons. First, the limitation 
of five registers allocated to the programmer by the PL/ATAC 
conpiler leads to a significant increase in program execution 
time as compared to a program written in ATAC assembly lancuage. 
As the programmer gains experience with the PL/ATAC languace, 
this limitation of the compiler may be overcome somewhat by 
programmer ‘tricks’ to gain use of registers. 

The lack of a double precision capability in the current 
vexsion of the PL/ATAC compiler installed on the IBM 360/57 
at the Naval Postgraduate School created a severe handicap 
in Writing the Pulse Train Separator Simulation Program and 
EMecOmCOMtGI nites tO program inefficiency. Each double 
precision calculation performed requires a program segment cf 
assembly language and the housekeeping involved results in 
increased program execution time. 

Even with these limitations it is felt that the simulation 


programs presented in this work provide the framework for ESM 
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software development and applications programming at the 


aval Postgraduate School. 


C. RECOMMENDATIONS 
Based upon the experience and results obtained from this 
Study, recommendations for follow-on study in the area of 
ESM software development are as follows: 
1. Continue study of the Pulse Train Separator Program 


and optimize the basic algorithm. 


2. Implement the resulting improved algorithm in ATAC 
assembly language. 


3. Use the Data Simulation Program and optimized hardware 
Simulation with APSS to develop and optimize ESM signal 
processing algorithms. 


4. Continve use of the PL/ATAC compiler to develop ané 
test hardware simulation and other programs prior to 
implementation in ATAC assembly language. 
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APPENDIX A 
INSTRUCTIONS FOR INSTALLATION AND USE OF THE ATAC 


PROGeoironG sUPPORT SYSTEM (APSS) ON THE IBM 360/67 
OPERATING SYSTEM AY THE NAVAL POSTGRADUATE SCHOOL 


Pv. INSTALLATION 
The ATAC Programming Support System is currently loaded 
Saecne 2321 Data Cell. Data Set names and their locations 


ace as fLoLlLows: 


ipo LOLS. RDS.LOADLIS ee ae O07 
Zee 5.616. RDS.APSS ~) CEL 005 
Se SLOlS.RDS.PLATAC = (00S) OEM 


In the event of failure of the Data Cell these data sets may 
be reloaded from the APSS source tape stored at the W.R. 
Church Computer Center under the name ‘ATI.006'. 

PibeomilmGiethiSswtape Contains a description of the tape. 
File #2 contains examples of the JCL required tc create and 
run APSS. The insStructions contained in this file are super- 
seded by the instructions contained here. The following JC 
is required to list these files: 


//(STANDARD JOB CARD) 

//EXEC PGM-IEBGENER 

//SYSPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT2 DD SYSOUT=A, DCPR=(LRECL=80 ,RECFM=F , BLKSIZE=80) 
LABEL (1,SL, ,1N),DCB= (BLXSIZE=800 ,LRECL=80, 
PECEM-F Sa i VOL—SCEhR-ATIO06 ,DSN=FILEL 

// EXEC PGM=IEBGENER 

J/BSISPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT2 DD SYSOUT=A , DCB= (LRECL=$89 , RECFM=F ,BLKSIZE=80) 

aoe Lemme UN td tT—(240C> DEFER) ,DISP=(OLD,PASS) , 
Pee ioe) peB=(BLKS1Ze2=800 ,LRECL=80) , 
RECFM=FB , VOL=SER=AT1I006 ,DSN=FILE2 


CW 
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File #3 of the source tape contains the object module 
ieoraryeot Aros and 1S requibed to build the APSS Instruction 
Level Simulator. File #4 contains the PL/ATAC Compiler object 
module. The following steps are required to build APSS from 
files 3 and 4. 

lL. Read files 3 and 4 from tape to disc storage. Special 
arrangements should be made to insure that the disc pack used 
has sufficient unallocated space for this step and the 
following step of this procedure. SPOOL3 has been made 
available in the past and is the disc used in the JCL of this 
example. The load modules from file 3 are moved into a disc 
file named 'S1614.RDS.LOADLIB'. File 4 1S moved into a Gisc 
fesiamnea “Sl6l5.RDS.PLATAC'. The utility control statements 
may be changed to rename the files as they are being loaded 
onto the disc. The JCL statements required are as follows: 


//(STANDARD JOB CARD) 

//STEP1 EXEC PGM=IEIIMOVE, REGION=1LOOK 

//SYSPRINT DD SYSOUT=A 

Pe evcuUrtl DD UNET=SYSDA,SPACE=(TRK, (40) , ,CONTIG) 

V7 ie DD) UNITT=2314 ,DSN=S1615.RDS.LOADLIB, 

Ve PACE =(URK, (50,10,10),,CONTIG), 

//DISP= (NEW , KEEP) , VOL=SER=SPOCOL3 

J/TITAPE DD UNIT=(2400,,DEFER) ,DISP=(NEW,PASS), 

fepeh=(s-Su,, IN) , 

V7 pee bEN=2,BiLKSILZE—800,LRECL=80 ,RECFM=FB), 
VOL=SER=ATI006 

//SYSIN DD * 

C6rr Pps -A7)  APSS. LOADLIB, £TO=2 314=SPOOL3 ,TODD=DA1, 
FROMDD=T1TAPE, * 

FROM=2 400=(ATIOCOG6,3) , RENAME=S1615.RDS.LOADLIB 

COPY DSNAME=ATI .APSS.PLATAC ,TO=2314=SPOOL3, 
FROMDD=TITAPE , TODD=DA2, * 

FROM=2400=(ATIO0O6,4) , RENAME=S1615.RDS.PLATAC 

77a? DD UNTT=2 314, pSN=S1615.RDS.PLATAC ,SPACE=(TRK, 
C20) REE); 

//DISP= (NEW,KEEP) , VOL=SER=SPOOL3 , LABEL=RETPD=360 

Ve * 
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Zee Une@sidbrary (named “S1615.RDS.LOADLIB') now on the 
disc must be built into an absolute load module of the entire 
APSS system. The module thus built (named 'S1615.RDS.APSS') 
will be stored permanently on cell 005 of the data cell. This 
file may also be renamed as desired. The following JCL state- 
ments are required. 


//{STANDARD JOB CARD) 

77/BULLD EXEC PGM=[EWL,REGION=L50K, 

Py Pak = Oviny, XREP, WET ,LIST,SIZE=(256K,20480) ' 

//SYSPRINT DD SYSOUT=A 

//LIBRARY DD DSN=S1615.RDS.LOADLIB,UNIT=2314, 
VOL=SER=SPOOL3 , DISP=SHR 

//SYSLIB DD DSNAME=3YS1.FORTLIB,DISP=SHR 

//SYSLMOD DD DSNAME=S1615.RDS.APSS,UNIT=2321, 
Ol=SER-CELOOS , 

//DISP= (NEW , KEEP) , LABEL=RETPD=360, 

Vy OPACE=(TRK, (300,20,2) ,RLSE) 

P7oYoUT) DD UNIT—SYSDA,SPACE=(TRK, (19,19) ,,CONTIG), 
SEP=SYSLMOD 

//SYSLIM DD* 

INCLUDE LIBRARYEPMUDLN 

CHANGE MSIMEIHESAPDN 

INCLUDE LIBRARY (APSSMON) 

INCLUDE LIBRARY (MSIM4,ASEM5,SIM16A,SIMTRi ,SIMIO1,GUL) 

INCLUDE LIBRARYEXPLMONN 

OVERLAY Al 

INSERT MSIMUL, *MSIMULA, IHENTRY , IHESAP 

INSERT MINT, IN, OUT 

INSERT IHEDBN, IHEXID 

INSERT IHEBSM, IHECSM 

INSERT IHEBSK, IHEIOX, IHEIOP, IHEDID, IHEDOB 

INSERT IHEBIB, IHEDOA, IHEIOB 

INSERT IHEIOA, IHEOCL 

INSERT IHEBSD, IHEBSF 

INSERT IHEJXS 

iio eli eOspe IHROST, ITHEBSI 

INSERT IHEVPE, IHEDMA, IHEVFB 

INSERT IHEDNC, IHEVFD, IHEVFA, IHEVPD, IHEVPB, IHEVSC 

INSERT IHEVSD, IHEVFE, IHEDCN, IHEUPB 

INSERT IHEVFC, IHEVPE, IHEVPG, IHEVQB, IHEVQC 

INSERTIHEABN, IHEIOD, IHEIOF, IHEPRT, IHEVQA, IHESPRT 

INSERT IHEBEG, IHEERR, IHESIZ. 

INSERT MISET 

OVERLAY Al 

INSERT ASSEM, REWIND, REW72, DSKOUT, CARDIN, DISKIN, 

ERPRT, PRIADD 


oe. 





INSERT WRDATA, PRICOM, PRINOR, WRITEX, REFTIT, PREF, 
Bie Ts 

OVERLAY Al 

INSERT PARMRD, PRESIM 

OVERLAY Al 

INSERT SMLTR, NRMTRM, STRTSM, TPAGE 

INSERT RDCRD, ABNRMT, ARTHFP, TRACE, HGRAM, HGRAMI, 
HGRAMS 

INSERT IOINIT, ACTIVE, SIMTIM 

INSERT DEVDTA, ACT, TIME, INT, RAND, DEADT, DEBUG 

OVERLAY A2 

INSERT LEVEL, DMAIOI, DMAIOA, DMAIOD, RIOIO 

TNGERT REMACT, DMA, DMATM, RIO, RIOTM, RIOINT, INTOLY, 
DMAIN? 

INSERT DTRAN, PUTACT, RANDOM 

OVERLAY Al 

INSERT HGPRNT 

OVERLAY Al 

INSERT LAODER 

OVERLAY Al 

GNSERT LINK, ENTEXT, SLLH 

OVERLAY Al 

INSERT PLATAC, IOPACK 

OVERLAY SOBJECTEREGIONN 

imicoeRT OBJECT, INIT, LIB, RCALPH, RCHEX, RCINT 

INSERT MDATE 

OVERLAY SDUMPEREGIONN 

INSERT SMDUMP, PAGE 

ENTRY MAIN 

NAME APSS 

/* 


3. The PL/ATAC compiler object module on the disc must 
be moved to the data cell. 


//(STANDARD JOB CARD) 

//EXEC PGM=IEBGENER 

//SYSPRINT DD SYSOUT=A 

PoerouT ll DD DISP=SHR,UNIT=2314,VOL=SBR=SPOOL3, 
//DSN=S1615.RDS.PLATAC 

//SYSUT2 DD DISP=(NEW,KEEP) ,UNIT=2321,VOL=SER=CELOOL, 
//DCB= (RECFM=FT ,BLKSIZE=7188) ,LABEL=RETPD=360, 

Poe er (7lse.(20,0)) ,DSGN=S1615.RDS.PLATAC 

J//SYSIN DD *# 

GENERATE 


/* 
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(leo m into maGy ool SakbDo,.LOADLIB') must now be moved 
from the disc to the data cell for permanent storage. This 
library contains modules necessary for successful execution 
of the PL/ATAC compiler. The required JCL statements to 
accomplish this step are as follows: 


//(STANDARD JOB CARD) 

//EXEC PGM=IEBCOPY 

//SYSPRINT DD SYSOUT=A 

//SYSUT1 DD DISP=SHR,UNIT=2314,VOL=SER=SPOOL3, 
V7 OSN=S1615.RDS.LOADLIB 

//SYSUT2 DD DISP=(NEW,KEEP) , UNIT=2321,VOL=SER=CELO07, 
//DSW=S1615.RDS.LOADLIB, SPACE=(7294,(61,0,14)), 
/ /DCB= (RECFM=UT , BLKSIZE=7294) , LABEL=RETPD=360 
Py eYSUT3 DD UNIT=SYSDA ,SPACE=(TRK, (20,5) ) 
Pyeyvour4a DD UNIT=SYSDA,SPACE=(TRK, (20,5)) 
//SYSIN DD * 

GOPY OUTDD=SYSUT2,IWDD=SYSUTI1 


/* 


Eee USER INSTRUCTION 
The JCL required to compile and execute PI/ATAC programs 
is as follows: 


//(STANDARD JCB CARD) 

//PLATAC EXEC PGM=XPLMON , PARM='ALTER' , REGION=200K 

//STEPLIB DD DSN=S1615.RDS.LOADLIB, UNIT=2321,D1ISP=SHR, 
VOL=SER=CEL007 

//PROGRAM DD DSN=S1615.RDS.PLATAC, UNIT=2 321 ,VOL=SER=CELOOL, 
DISP=SHR 

-//SYSPRINT DD SYSOUT=A 

//FTO5F001 DD DDNAME=SYSIN 

V7ELOGEOOL DD SYSOUT=A 

J/ET2Z7F001 DD DSN=&&PLCODE,SPACE=(TRK, (20,10) ,RLSE), 

y/o DCB= (BLKSIZE=3520, LRECL=80, RECFM=FB) , DISP=(NEW,PASS,DELETE) , 

Ve UNIT=SYSDA 

//SYSIN DD * 


PL/ATAC Program 





//APSS EXEC PGM=APSS REGION=200K 
//STEPLIB DD DSN=S1615.RDS.APSS ,UNIT=2321,DISP=(OLD,KEEP) , 
//  VOL=SER=CEL005 , DCB=BLKSIZE=400 


//FTO5F001 
//FTO6FO01 
//FT08FO01 
//FTO9FOOL 


P77 ETAL O01 


V/ET20F.001 


DD 


DD 


DDNAME=SYSIN 

SYSOUT=A , DCB= (RECFM=FBA , LRECL=133, BLKSIZE=3325) 
UNIT=SYSDA , SPACE= (CYL, 1) , DCB=BUFNO=1 
UNIT=SYSDA ,SPACE= (CYL, (7,2)), 

DCB= (RECFM=VBS . BLKSIZE=7180 , LRECL=92 , BUFNO=1) 
UNIT=—¢ ~SmA SPAGE= (CYL, (7,2)), 

DCB= (RECFM=VBS , BLKSIZE=4204 , LRECL=42 , BUFNO=1) 
Pim ov cuAmeorpncCe =(GVL. (7,2) )— 


DCB= (BUFNO=1, RECFM=VBS , BLKSIZE=2004 , LRECL=500) 

//FT27F001 DD DSN=&&PLCODE, UNIT=SYSDA, DISP= (OLD, DELETE) , 
DCB=BUFNO=1 

//SYSPRINT DD SYSOUT=A , DCB= (RECFM=FBA , LRECL=133 , BLKSIZE= 3325) 


Y/SYXSIN DD 
SEND 
Vas 


* 
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APPENDIX B 


DATA SIMULATION PROGRAM FORTRAN LISTING AND OUTPUT 


A description of program variables is included in the 


Data Simulation Program listing. The input data deck setup 


is also included in the listing. 
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PULSE TRAIN SEPARATOR SIMULATION PROGRAM PL/ATAC LISTING 


Program variables are deScribed as they are declared in 
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